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Much of atmospheric and oceanic transport is associated with coherent struc- 
tures. Lagrangian methods are emerging as optimal tools for their identification and 
analysis. An important Lagrangian technique which is starting to be widely used 
in oceanography is that of Finite-Size Lyapunov Exponents (FSLEs). Despite this 
growing relevance there are still many open questions concerning the reliability of 
the FSLEs in order to analyse the ocean dynamics. In particular, it is still unclear 
how robust they are when confronted with real data. In this paper we analyze the 
effect on this Lagrangian technique of the two most important effects when facing 
real data, namely noise and dynamics of unsolved scales. Our results, using as a 
benchmarch data from a primitive numerical model of the Mediterranean Sea, show 
that even when some dynamics is missed the FSLEs results still give an accurate 
picture of the oceanic transport properties. 
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Lagrangian viewpoint [5l [H [THl [22l [24] . Lagrangian diagnostics exploit the spatio- 
temporal variability of the velocity field by following fiuid particle trajectories, in con- 
trast with Eulerian diagnostics, which analyze only frozen snapshots of data. Among 
Lagrangian techniques the most used ones involve the computation of local Lyapunov 
exponents (LLE) which measure the relative dispersion of transported particles [H EJ HI 
[T3l [m [IBJ. LLEs give information on dispersion processes but also, and even more im- 
portantly, can be used to detect and visualize Lagrangian Coherent Structures (LCSs) in 
the fiow like vortices, barriers to transport, fronts, etc. [3l [T5l fT9l - [2n [27] . 

The standard definition of Lyapunov exponents [26j involves a double limit, in which 
infinitely-long times and infinitesimal initial separations are taken. These limits can not 
be practically implemented when dealing with realistic fiows of geophysical origin. Over 
real data, LLEs are defined by relaxing some of the limit procedures. In finite-time 
Lyapunov exponents (FTLE) [211 ES] trajectory separations are computed starting still 
at infinitesimal distance, but only for a finite time. In the case of finite-size Lyapunov 
exponents (FSLE), the key tool used in our work [21 [4], one computes the time which is 
taken for two trajectories, initially separated by a finite distance, to reach a larger final 
finite distance. 

FSLEs are attracting the attention of the oceanographic community [H [fij [Tj, [HI [EJ [HI 
\2M [281 [30 j. The main reasons for this interest, within the framework of their ability for 
studying LCSs and dispersion processes, are the following: a) they identify and display 
the dynamical structures in the fiow that strongly organizes fiuid motion (the above 
mentioned LCSs, which are defined as ridges of the FSLEs fields); b) they are relatively 
easy to compute; c) they provide extra information on characteristics time-scales for the 
dynamics; and d) they are able to reveal oceanic structures below the nominal resolution 
of the velocity field being analyzed. In addition, FSLEs are specially suited to analyze 
transport in closed areas [Ij. 

Despite the growing number of applications of FSLEs, a rigorous analysis of many of 
their properties is still lacking. There are two main concerns before applying FSLE to real 
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data, namely the effect of noise and the role of observation scale. Concerning noise, real 
data are discretized and noisy, and this can affect the reliability of FSLE-based diagnostics 
(recent related studies for FTLE can be consulted in [27]). Concerning scale properties, 
FSLEs can be obtained over a grid finer than that of the data. This enables to study 
submesoscale processes under the typical mesoscales (below 10 kilometers) that nowadays 
provide, for example, altimetry data [7l[3^. But then the question is if finer-grid LCSs are 
meaningful or just an artifact. On the other side, we can have access to a limited-resolution 
velocity field, and then ask ourselves if any refinement in the velocity grid (by improved 
data acquisition, for instance) is going to modify our previous assessment of LCS at the 
rougher scale. The main objective of this paper is to address these questions, in particular 
with reference to their potential applications into ocean dynamics. A related study of the 
sensitivity of relative dispersion of particles statistics, when the spatial resolution of the 
velocity field changes, has appeared recently [28] . 



The benchmark for the study of FSLEs properties used in this work is a the two- 
dimensional velocity field of the marine surface obtained from a numerical model of the 
Mediterranean Sea. The first half of the paper is devoted to study what we have just 
called scale properties of the FSLE field. By changing gradually the resolution of the grid 
on which they are computed, we will show that they have typical multifractal properties. 
This means, in particular, that FSLEs obtained for a finer resolution than that provided 
by the data provides non-artificial information. Subsequently, we will consider a somewhat 
opposite case, i.e., what happens to the FSLEs if the velocity-data grid is changed. Their 
robustness under data-resolution transformations will be discussed. The second half of 
the work will analyze the effect of noise. Again, two different scenarios are considered: 
a) uncertainties in the velocity data, and b) noise in the particle trajectories. FSLEs will 
be shown to be robust against these two sources of error, and the reasons for this will be 
discussed. 
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I. DESCRIPTION OF THE DATA 

We analyze a velocity dataset generated with the DieCast (Dietrich for Center Air Sea 
Technology) numerical ocean model adapted to the Mediterranean Sea [9j. The dataset 
has been already used in previous Lagrangian studies [SJ [23l [29j . DieCast is a primitive- 
equation, z-level, finite diflFerence ocean model which uses the hydrostatic, incompressible 
and rigid lid approximations. At each grid point, horizontal resolution is the same in 
both the longitudinal, 0, and latitudinal. A, directions, with resolutions A(/) =1/8^ and 
AA = A(/)cos A. Vertical resolution is variable with 30 layers. Annual climatologic forcing 
is used, so that it is enough to keep a temporal resolution of one day. We will use velocity 
data corresponding to the second layer, which has its center at a depth of 16 meters. This 
sub-surface layer is representative of the marine surface circulation and is not directly 
driven by wind. We have recorded daily velocities for five years, and concentrate our 
work in the area of the Balearic Sea. In Figure [T] we show a snapshot of the velocity field 
from the DieCAST model. 
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FIG. 1. Snapshot of the surface velocity field of the Balearic Sea corresponding to day 640 in 
the DieCAST simulation. 
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II. DEFINITION AND IMPLEMENTATION OF FSLES 

FSLEs provide a measure of dispersion as a function of the spatial resolution, serving 
to isolate the different regimes corresponding to different length scales of oceanic flows, 
as well as identifying the LCSs present in the data. To calculate the FSLEs we have to 
know the trajectories of fluid particles, which are computed by integrating the equations 
of motion for which we need the velocity data (u^v). FSLE are computed from r, the 
time required for two particles of fluid (one of them placed at {x^y)) to separate from an 
initial (at time t) distance of 5o to a final distance of , as: 

1 5f 

A{x,y,t,5o,Sf) = -log/. (1) 

T Oq 

In principle obtaining A(x, t, So^6f) would imply to consider all the trajectories start- 
ing from points at distance Sq from our basis point; in practice, when confronted with 
regular, discretized grids, only the four closest neighbors are considered. It is convenient 
to choose Sq to be the intergrid spacing among the points on which the FSLEs will be 
computed, i.e., it is the resolution of the "FSLE grid". The details of the calculation of 
the FSLEs are in Appendix |A| 

III. EFFECT OF SAMPLING SCALE ON FSLES 

Notice that we distinguish two types of grids, namely the FSLE grid (in the following 
"F-grid") and that where the velocity field is given (called velocity grid or "V-grid"). 
These two grids need not to coincide. We explore first the effect of changing the resolution 
of FSLE grid. 

In Figure [2] we show an example of the FSLEs derived using four different F-grids with 
increasing resolutions. Visually, as the resolution of the F-grid is increased the structures 
already observed in the coarser version are kept, just increasing their detail level. In 
addition, when resolution is increased, new less active structures are appearing in areas 
previously regarded as almost inactive. Taking a F-grid finer than the associated V-grid 
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would make no sense if FSLE was an Eulerian measure obtained from single snapshots 
of the velocity field. But FSLE is a Lagrangian measure, i.e. they are computed using 
trajectories which integrate information on the history of the velocity field. This allows 
capturing the effects of the large scales on scales smaller than the V-grid. This does not 
mean that we reconstruct all the eflFects taking place at the smaller scales, but only the 
ones that have been originated by the relatively large velocity features which are resolved 
on the V-grid. 




0,00 0,10 0.20 0.30 O.+O 0,50 0,B0 



FIG. 2. Snapshots of FSLEs backward in time starting from day 640 at different F-grid resolu- 
tions: a) 60 = 1/8°, b) 5o = 1/16°, c) So = 1/32°, d) So = 1/64°. In all of them we take Sf = 1°. 
The values in the color bar have units of day~^. 

The main structures of the fiow, which are essentially filaments, become finer as reso- 
lution is increased, behaving much like the geometrical persistence of a fractal interface 
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[8]. The question naturally rises about the possible multifractal character of FSLEs. 
Multifractality is a property characteristic to turbulent flows, and it is associated to the 
development of a complex hierarchy accross which energy is transmitted and utterly dis- 
sipated [TOl [33]. The study of the scaling properties of the distribution of FSLEs at the 
different resolution scales reveal that they are multifractal (see Appendix [B]). This im- 
plies that changes in scale are accounted for by a well-deflned transformation, namely a 
cascade multiplicative process ^01 |33j. It also implies that information is hierarchized 
[3T J and so what is obtained flrst, at the coarsest scales, is the most relevant information. 
Due to multifractality, small-scale structures, as unveiled by the FSLEs, with typical sizes 
smaller than that of the velocity resolution, are determined by the larger ones and the 
multi-scale invariance properties. Thus, no artiflciality is induced by this calculation and 
the robustness of FSLE analysis under changes in scale is conflrmed. 

A diflFerent question concerns the robustness of FSLEs when the V-grid is changed. 
For instance, when a diagnosis is obtained with a low-resolution velocity fleld, is this 
diagnosis compatible with a later improved observation of the velocity? The answer 
is yes. In Figure |3] we show the FSLEs obtained at a flxed F-grid resolution of 1/8° for 
varying velocity resolutions (namely, 1/8°, 1/4° and 1/2°). The change of resolution of the 
velocity grid is performed as indicated in Appendix [Cl We observe that the global features 
observed with the coarser resolution V-grid are kept when this is reflned. Obviously, as 
the velocity fleld is described with enhanced resolution new details (with short-range effect 
on the flow structure and so no contradicting the large-scale picture) become apparent in 
the FSLEs. However, the effect of reflning the V-grid is not only introducing new small- 
scale structures: there is a consistent increase in the values of the FSLEs as the V-grid is 
reflned. The histograms of the FSLEs, Aj, for a given velocity resolution conditioned by 
the FSLEs, Ac, obtained from a coarser velocity fleld are shown in Figure |4} The modal 
line (the line of maximum conditioned probability) is close to a straight line. The best 
linear regression flts are Ai/8° = I.O8A1/40 +0.05 (correlation coefflcient, p = 0.69) and 
Ai/40 = 0.99Ai/2° + 0.04 (p = 0.69). According to these results, we can approximate the 
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finer FSLEs Aj in terms of tlie coarser FSLEs Ac as: 

A;(f) = Ac(f) + AAjc, (2) 

wliere the quantity AAjc determines the contribution to the FSLE by the small-scale 
variations in velocity not accounted for by the lower resolution version of this field. It is 
hence independent of A^ so the intercept of the vertical axis with the linear regression 
equals the mean of this quantity. 

A linear dependence of A/ with Ac when scale is changed implies that FSLEs follow a 
multiplicative cascade [321 an essential ingredient in multifractal systems which gives 
further confirmation to our previous results. We have hence shown that a) the dependence 
of FSLEs on both types of scale parameters reveals a multifractal structure; b) what is 
diagnosed at the coarser scales is still valid when scale is refined (although as V-grid 
resolution is increased the reference level of FSLEs is increased by a constant). 



IV. EFFECT OF NOISE ON FSLES 



We compute the FSLEs after applying a random perturbation to all components of the 
velocity field. The velocity is changed from (u^v) to {u\v')^ with u'{x^t) = 2i(x, + 
ar]x{^^t)) and v'{^^t) = '^(x, + ar]y{-x^t)). {/^^^(x, t), 7y^(x, t)} are sets of Gaussian 
random numbers of zero mean and unit variance, a measures the relative size of the 
perturbation (it gives the ratio of the mean amplitude of noise with respect to mean 
amplitude of the velocity). We introduce three different kinds of error: uncorrelated noise, 
i.e. different and uncorrelated values of {/^^^(x, t), /^^/(x, t)} for each x and t; correlated in 
time and uncorrelated in space (uncorrelated for different x but the same values at given 
X for different t); and correlated in space and uncorrelated in time (uncorrelated values 
for different t, but the same values for different x at fixed t). Note that the perturbation is 
proportional to the original velocity. Figure [5] shows a snapshot of FSLEs for the velocity 
field perturbed by uncorrelated noise of a relative size a = 10, i.e., noise is 10 times larger 
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FIG. 3. Snapshots of FSLEs backward in time starting from day 640 at different initial V-grid 
resolutions: a)Ao = 1/8°, b) Aq = 1/4°, c) Aq = 1/2°. In all of them we take the same F-grid 
resolution of 5o = 1/8°, and 5j = 1°. The color bar has units of day"-*^. 

than the amplitude of the initial velocity field. The computed Lagrangian structures look 
rather the same, despite the large size of the perturbation introduced. 

In order to quantify the influence of the velocity perturbation in the FSLE calculation 
we compute the relative error (RE) of perturbed FSLEs with respect to unperturbed 
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FIG. 4. Probability distributions (coded as grey levels) of FSLEs derived at a coarse velocity 
resolution (A^ vertical axis) conditioned by a finer velocity grid (Aj, horizontal axis). The range 
of values of both axes run linearly from to 0.5 day~^. To obtain a statistics large enough, 
we have considered 30 FSLEs snapshots, starting form t = 640 up to t = 1075 days, in steps 
of 15 days. The brightest color (pure white) corresponds to the maximum probability at each 
column; the darkest color (pure black) corresponds to zero. Left panel: FLSEs derived from 
1/4° velocities conditioned by FSLEs at original 1/8° velocities. Right panel: FSLEs from 1/2° 
velocities conditioned by FSLEs from 1/4° velocities. 

FSLEs, < e{t) >, at a given instant of time, and then averaging in time (we have M = 100 
snapshots : t = ti, ...^tM) as: 

A(x, ti) and A'^(x, ti) are the FSLEs fields without and with inclusion of the perturbation 
in the velocity data, respectively. The sum over x runs over the = 2679 spatial points. 
Figure [6] shows the average RE as a function oi a. It must be remarked that the RE has 
always small values: even for a = 10 the RE remains smaller than 0.23 for the three kinds 
of noise. To get an idea of how relevant these quantities are, we have computed the RE 
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FIG. 5. Snapshots of FSLEs calculated backwards in time starting from day 600 at fixed spatial 
resolution (Sq = 1/64°), and at different a: a)a = 0, b) a = 10. In both of them we take 6f = 1°. 
The color bar has units of day~^. Initial conditions for which the separation 6f has not been 
reached after 600 days are assigned a value A = 0. 



of shuffled FSLEs (permuting locations at random) with respect to the original ones, and 
obtained a value of 1.143. FSLEs are thus robust against uncorrelated noise; the reason 
is the averaging effect produced when computing them by integrating over trajectories 
which extend in time and space, that tends to cancellate random, uncorrelated errors. 

Now we proceed by adding noise to the particle trajectories. This is a simplified 
way of including unresolved small scales in the Lagrangian computations [12]. To be 



precise we solve numerically the system of Equations (Dl) and (D2) (see Appendix D), 
where a stochastic term with a Gaussian random number and an effective eddy-diffusion, 
has been added. For the diffusivity we use Okubo's empirical formula [25j, which 
relates the effective eddy-diffusion, D in rri^ / with the spatial scale, / in meters: D{1) = 
2.055 10~^ l^'^^. If we take / = 12 /cm, which is the approximate length corresponding to 
the 1/8° DieGAST resolution at Mediterranean latitudes, we obtain ^ 10 m^s~^ = L^o- 

In Figure [T] we show particle trajectories without (top panel) and with (bottom panel) 
eddy diffusion. As expected diffusion introduces small scale irregularities on the trajecto- 
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FIG. 6. Relative error < e{t) > of the FSLE fields for different perturbation intensity a in the 
velocity data. Solid line is for uncorrelated noise in space and time, dashed-dotted line is for 
uncorrelated noise in time and correlated in space, and dotted line is for uncorrelated noise in 
space and correlated in time. < e(t) > is obtained by averaging the RE in 100 snapshots (see 
Eq. ([3])). The error bar is the statistical error of the temporal average < e{t) >. Left: spatial 
resolution Sq = 1/8°. Right: spatial resolution Sq = 1/64°. In all calculations we take Sf = 1° 

ries, and also substantial dispersion at large scales. In Figure |8| FSLEs with Sq = 1/64° 
D = m^s~^ and D = 0.9 m^s~^ (obtained for this scale by Okubo's formula for spatial 
resolution) are shown. We can see that the main mesoscale structures are maintained, but 
small-scale filament al structures are lost since filaments become blurred. This is some- 
how expected a priory because diffusion introduces a new length scale Id proportional to 
\fD. A pointwise comparison of noiseless and noise- affected FSLEs makes no sense, since 
the noise-induced blurring disperses FSLEs values specially at places with low values, see 
Figure [8| But for Lagrangian diagnostics high FSLE values are much more relevant, so 
we compute the error restricted to the places where A > 0.2 for D = 0. Left panel of 
Figure [9] shows the RE respect to the = case (applying Equation ([3])) for different 
values of D. The RE monotonously increases with but remains smaller than 0.6 for 
the largest value of D considered. Again, to get an idea of how relevant these relative 
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errors are, we have to compare them with the RE of shuffled FSLEs, of value 1.143. 

As a matter of fact, diffusion introduces an effective observation scale, and one should 
not go beyond that limit to obtain senseful results; this is illustrated in the right panel 
of Figure [9| As shown in the figure, for fixed eddy diffusivity (in the case of the figure. 
Do = 10m^5~^), when 60 becomes greater the error diminishes. Hence, a fixed diffusion 
will eventually be negligible at a scale large enough, in this way determining an observation 
scale. A completely different situation is given when the eddy-diffusion depends on the 
scale according to Okubo's formula (in our case, at Sq = 1/16° D = 4, 5m^5~-^; at 5o = 
1/32° D = 2m^s-^] and at So = 1/64°, D = 0.9m^s-^). Now < e{t) > takes a constant 
value close to 0.45, meaning that Okubo's diffusion behaves the same at all scales. This 
is expected since Okubo's law is based on the hypothesis that unsolved scales act as 
diffusers, like in turbulence [lOj. Our result is consistent with the ideas behind Okubo's 
hypothesis. 

We can characterize the effective diffusion scale from the properties of the histograms. 



In Figure [10] we present the histograms of FSLEs at different F-grid resolutions, and 
including or not eddy diffusion which takes always the same value Dq = 10m^s~^. For 
(^0 = 1/8° the histograms with and without diffusion are almost coincident. This is due to 
the fact that the value of diffusion we are using is the one corresponding, by the Okubo's 
formula, to 1/8°. i.e., we are parameterizing turbulence below 1/8°, and this has no 
effects on the FSLE computations if the minimum scale considered is also 1/8. However, 
this behavior is different for smaller 60 (always keeping Dq = 10rn^s~^). The histograms 
for 5o = 1/16,1/64°, with and without diffusion, are clearly different, those including 
diffusion becoming closer to the histogram for Sq = 1/8°. 



V. CONCLUSIONS 

In this paper, we have analyzed the sensibility of FSLE-based analyses for the diagnos- 
tic of Lagrangian properties of the ocean (most notably, horizontal mixing and dispersion). 
Our sensibility tests include the two most important effects when facing real data, namely 
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dynamics of unsolved scales and of noise. Our results show that even if some dynamics 
are missed (because of lack of sampling or inaccuracy of any kind in the measurements) 
FLSEs results would still give an accurate picture of Lagrangian properties, valid for the 
solved scales. This does not mean that scale and/or noise leave FSLEs unaffected, but 
the way in which they modify this Lagrangian diagnostics can be properly accounted. 

Lagrangian methods provide answers to problems which have a deep impact on risk 
management (e.g. control of pollutant dispersion) as well as on ecosystem analysis (e.g. 
tracking nutrient mixing and transport, identifying the role of horizontal mixing in pri- 
mary productivity). They utterly will give hints about energy exchanges in the upper 
ocean and will help in understanding processes driving global change in the oceans. The 
use of Lagrangian techniques for the assessment of the transport and mixing properties 
of the ocean has grown in importance in the latest years, with increasing efforts devoted 
to the implementation of appropriate techniques but few studies on the validity of the 
results when real data, affected by realistic constraints, are used. Our work will serve 
to unify and interpret the analyses provided by Lagrangian methods when real data are 
processed. 
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Appendix A: Calculation of FSLE. 

It is natural to choose the initial points {x^y) on the nodes of a grid with lattice 
spacing coincident with the initial separation of fluid particles 5o. In this way one obtains 
maps of values of A at a spatial resolution that will coincide with Sq. To compute A we 
need to know the trajectories of the particles. The equations of motion that describe the 
horizontal evolution of particle trajectories in the velocity field are: 
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dt ~ RcosX ' ^ ' 

dX v{(f>,X,t) 



dt R 



(A2) 



where u and v represent the zonal and meridional components of the surface velocity 
field coming from the simulations described above, R is the radius of the Earth (6400 km 



in our computations). Numerically we proceed integrating Eqs. (Al) and ( |A2| ) using a 
standard fourth-order Runge-Kutta scheme with an integration time step dt = 6 hours. 
Spatiotemporal interpolation of the velocity data is achieved by bilinear interpolation. 
Since this technique requires an equally spaced grid and this is not the case for the 
spherical coordinates (0, A), for which the grid is not uniformly spaced in the latitude 
coordinate, we first transform it to a new system (0, /x) where the grid turns out to be 
uniformly sampled [23] . The latitude A is related to the new coordinate fi according to: 

/i = log I sec A + tanA| (A3) 
Using the new coordinate variables the equations of motion become: 



dcj) 
~dJt 



Rco8 A(/x) 
Rco8 A(/x) ' 



(A4) 
(A5) 



and one can convert the fi values back to latitude A by inverting Eq. (A3): A = 7r/2 — 



2arctane ^. Once we integrate the equations of motion, Eqs. (A4) and (A5), we com 



pute the FSLEs applying Eq. ([T]) for the points of a lattice with spacing 5q. We will 
only compute FSLEs integrating backwards-in-time the particle trajectories, since LCSs 
associated to this has a direct physical interpretation, but all our results are similar for 
forward-in-time dynamics. 
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Appendix B: Multifractal character of FSLE. 

To study if FSLEs behave like multifractal systems we have computed their probability 
density distribution, P{6o^A)^ for different resolutions 60. It must follow that for any 
resolution scale Sq of FSLE grid we should observe: 



P(,5o,A)oc5o'"''^^\ (Bl) 

Notice that as d—D{A) is a positive quantity, P{So^ A) becomes smaller as 60 is reduced. 
In fact, a characteristic signature of multifractal scaling is having a scale-dependent his- 
togram which becomes more strongly peaked as the resolution scale becomes smaller. 
We know that there is a FSLE for each domain point, so we can normalize the FSLE 
distribution by its maximum (that will be attained for a value Ac) and then retrieve the 
associated singularity spectrum according to the following expression: 



1 P{So,A) 

D{A) = rf- ^^^^t^-) . (B2) 
log do 



In Figure top, we show the histograms (averaged over 30 snapshots distributed 
among 15 months) normalized to have the same unitary are. One can see that when the 
resolution gets finer the P{So^ A) narrows, and the peak height increases, in agreement with 



Equation (Bl). According to the definition of Microcanonical Multifractals [13, [33], the 



system will be multifractal if the curves D{A) estimated at diflFerent resolutions 5q using 



Equation (B2) are all equal. This can be observed in Figure 11, bottom. The collapse 
of the four curves is not perfect due to the lost of translational invariance produced by 
the small size of the domain -the Balearic basin- that we analyze. Thus, the interfaces 
of constant A values build an approximate multifractal hierarchy and generalized scale 
invariance is present in the FSLE field. 
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Appendix C: Change of resolution of the velocity grid. 

In order to reduce the resolution of a given velocity grid, the easiest way would be by 
subsampling the existing grid points or by block-averaging the values of the velocity and 
assign the result to the central grid point. However, in the case of a complex boundary 
such as the Mediterranean coast such a strategy is strongly inconvenient, as the coars- 
ening of the grid would imply to change land circulation barriers (islands, straits). The 
disappearance of a land barrier or the creation of a new one as a consequence of the 
coarsening would imply a dramatic change in the value of FSLEs at all points affected 
by the modified circulation; if the circulation patterns are rather complex almost every 
point could be affected. We have thus preferred to smooth the velocity with a convolution 
kernel weighted with a local normalization factor, and keeping the original resolution for 
the data so land barriers are equally well described than in the original data. 

We define the coarsening kernel of scale factor 5, Hh g , as I 

Ks{oo,y) = (CI) 

We disregard the introduction of a normalization factor at this point since we will need 
to normalize locally later. The coarsened version of the velocity vector would hence be 
given by the convolution of this kernel with the velocity, denoted by Kgi^v. A coarsening 
convolution kernel turns out to be convenient with almost horizontal fields, as the deriva- 
tives commute with the convolution operator so if \/v ^ hence His ^ ViT ~ 0. However, 
this coarsening scheme needs to be improved. By convention, we take the velocity v as 
zero over land points. For that reason, a simple convolution does not produce a correct 
coarsened version of the velocity because points close to land would experience a loss of 
energy. The easiest way to correct this is to normalize by the weight of the sea points. 
Let us first define the sea mask M{x,y) as 1 over the sea and over the land. The 
normalization weight is given by Kg ^ M. For points very far from the land, this weight is 
just the normalization of Kg. For points surrounded by land points the weight takes the 
contributions from sea points only. We thus define the coarsened version by a scale factor 
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s of the velocity, Vg^ as: 



Vg r- (C2) 



In Figure [12] we show two examples of coarsened velocities. We can see that typical 
circulation patterns are coarsened as s increases, while land obstacles are preserved. In 
fact, if Aq is the velocity resolution scale, the effective resolution scale of Vg is sAq (the 
nominal resolution, on the contrary, is the original one, Aq). 



Appendix D: Introducing noise in the particle's trajectories. 

In order to introduce noise in the particle's trajectories we resolve the following system 
of stochastic equations: 

dt i?cos(A) ^i?cos(A)^'^ ^' ^ ' 

dX _ v{ct),X,t) V2DUt) 

dt~ R ^ R ' ^^^^ 
^ = 1,2 are the components of a two-dimensional Gaussian white noise with zero 



mean and correlations < ^= ^ij^{t — t').Eqs. (Dl and D2) use a simple white 

noise added to the trajectories. A more realistic representation of small-scale Lagrangian 
dispersion in turbulent fields requires using other kinds of correlated noises [12j but, as we 
are interested in examining infiuences of the missing scales, it is convenient to use white 
noise, since this would represent the extreme case of very irregular trajectories which gives 
an upper bound to the effects of more realistic smoother small scales. Thus the tests 
presented here are similar to the ones considered before (perturbation of the velocity) 
when adding uncorrelated perturbations to the velocity, but here the perturbation acts 
at arbitrarily small scales, as appropriate for a turbulent field, instead of being smooth 
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below a cutoff scale, as appropriate for modelling observational errors. 
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FIG. 7. Trajectories of five particles without diffusion (top) and with diffusion (bottom). The 
difference in the initial positions of all five particles is about 0.06°, and we use these initial 
conditions in both computations. The trajectories were computed for 50 days of integration. 
We used the eddy-diffusion Dq ^ 10m?s~^ assigned by the Okubo formula to the resolution of 
the DieCAST model at Mediterranean latitudes. 
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FIG. 8. FSLEs computed backwards from day 500 at the same spatial resolution {60 = 1/64°), 
and for different eddy-diffusion values: a)i? = m?s~^ h) D — 0.9 m?s~^ . We take 5f — 1°. 
The color bar has units of day~^. Initial conditions for which the separation 6f has not been 
reached after 500 days are assigned a value A = 0. 
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FIG. 9. Left: Relative error < e(t) > of the FSLE at the different values of D in the particle 
trajectories, with respect to the D = case. Spatial resolution is = 1/8°? ^iid 5f = 1°. 

< e(t) > is obtained by temporally averaging the relative errors in 100 snapshots. The (small) 
error bars indicate the statistical error in the < e{t) > average. Right: Dotted line is the RE 

< e{t) > of the FSLE at different spatial resolution 5o and at one assigned eddy-diffusion D for 
every spatial resolution in the particle trajectories with respect to the D = case. Solid line is 
the RE < e{t) > of the FSLE at different spatial resolution 5o, and at the same eddy-diffusion 
Dq = Wm?s~^ in the particle trajectories with respect to the D = case. Dashed-dotted line is 
the RE of shuffled FSLE with respect to the original case {D = 0) at different spatial resolution. 

< e{t) > by temporally averaging the RE in 100 snapshots. The (small) error bar indicates the 
statistical error in the < e(t) > average. In all of them we take Sf = 1^. 
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FIG. 10. Comparison between probability density function for the FSLEs at different resolutions 
with different values of eddy- diffusion, and without diffusion. It is obtained from the temporal 
average (30 snapshots) of histograms. Solid line for = 1/8° with diffusion, dotted 5o = 1/8° 
without diffusion, dashed = 1/16° without diffusion, dashed-dotted line for = 1/64° with 
diffusion, and circle-line for = 1/64° without diffusion. 
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FIG. 11. Top: Comparison of the probability density functions P((5o,A) for the FSLEs at 
different resolutions. It is obtained from the temporal average (30 snapshots) of histograms. 
Dotted line is for Sq = 1/8°, dashed-dotted line Sq = 1/16°, dashed Sq = 1/32°, and solid line 
for ^0 = 1/64°. Bottom: D{A) for different values of ^o- Dotted for = 1/8°? dashed-dotted 
line So = 1/16°, dashed = 1/32°, and solid line line for = 1/64°. 
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FIG. 12. Top: Velocity field coarsened by a scale factor s = 2, for a equivalent resolution 
Sq = 1/4°. Bottom: Velocity field coarsened by a scale factor s = 4, for a equivalent resolution 
So = 1/2°. 



